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ABSTRACT 

We study the dynamics of the FLRW flat cosmological models in which the 
vacuum energy varies with time, A(t). In this model we find that the main 
cosmological functions such as the scale factor of the universe and the Hubble 
flow are defined in terms of exponential functions. Applying a joint likelihood 
analysis of the recent supernovae type la data, the Cosmic Microwave Back- 
ground shift parameter and the Baryonic Acoustic Oscillations traced by the 
Sloan Digital Sky Survey (SDSS) galaxies, we place tight constraints on the 
main cosmological parameters of the A(t) scenario. Also, we compare the A(t) 
model with the traditional A cosmology and we find that the former model 
provides a Hubble expansion which compares well with that of the A cos- 
mology. However, the A(t) scenario predicts stronger small scale dynamics, 
which implies a faster growth rate of perturbations with respect to the usual 
A-cosmology, despite the fact that they share the same equation of state pa- 
rameter. In this framework, we find that galaxy clusters in the A(t) model 
appear to form earlier than in the A model. 

Keywords: cosmology: theory - large-scale structure of universe 



1 INTRODUCTION 

The detailed analysis of the available high quality cos- 
mological observations (Riess et al. 1998; Perlmutter 
et al. 1999; Efstathiou et al. 2002; Basilakos & Plionis 
2005; Tegmark et al. 2006; Davis et al. 2007, Kowalski 
et al. 2008; Komatsu et al. 2009) have converged during 
the last decade towards a cosmic expansion history that 
involves a spatial flat geometry and a recent accelerat- 
ing expansion of the universe. This expansion has been 
attributed to an energy component (dark energy) with 
negative pressure which dominates the universe at late 
times and causes the observed accelerating expansion. 
The simplest type of dark energy corresponds to the 
cosmological constant (see for review Peebles & Ratra 
2003). The nature of the dark energy is still a mystery 
and indeed it is one of the most fundamental current 
problems in physics and cosmology. 

In the literature, there are many theoretical spec- 
ulations regarding the physics of the above exotic dark 
energy. The simplest approach is to consider a real 
scalar field <f> which rolls down the potential energy 
V(<j>) and therefore it could mimic the dark energy (Ra- 
tra & Peebles 1988; Weinberg 1989; Turner & White 
1997; Caldwell, Dave & Steinhardt 1998; Padmanabhan 
2003). Alternatively, Ozer & Taha (1987) proposed a 
different scenario in which a time varying A parame- 
ter could be a possible candidate for the dark energy 
(see also Bertolami 1986; Freese et al. 1987; Peebles & 



Ratra 1988; Carvalho, Lima & Waga 1992; Overduin & 
Cooperstock 1998; Bertolami & Martins 2000; Alcaniz 
& Maia 2003; Opher & Pellison 2004; Bauer 2005; Bar- 
row & Clifton 2006; Montenegro & Carneiro 2007 and 
references therein) . In this cosmological model the dark 
energy equation of state parameter w = Pde / Pde, is 
strictly equal to -1, but the vacuum energy density (or 
A) varies with time. It is interesting to mention here that 
the renormalization group (RG) in quantum field theory 
(Shapiro & Sola 2000; Babic et al. 2002) provides a time 
varying vacuum, in which the A component evolves as 
~ H 2 (t) [see Grande, Sola & Stefancic 2006] , where H 
is the Hubble parameter. On the other hand, based on 
the holographic principle (Bousso 2002; Padmanabhan 
2005) one can prove that A ~ H 4 . 

However, in the A(t) cosmological model there is a 
coupling between the time-dependent vacuum and mat- 
ter (Carneiro et al. 2008). In particular, using the com- 
bination of the conservation of the total energy with 
the variation of the vacuum energy one can prove that 
the A(t) model provides either a process of a particle 
production or the mass of the dark matter particles in- 
creases. The latter general properties can be explained 
within the framework of the interacting dark energy 
models (Alcaniz & Lima 2005 and references therein). 
We would like to stress here that most of the recent 
papers in dark energy studies are based on the assump- 
tion that the dark energy evolves independently of the 
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dark matter. Of course, the unknown nature of both 
dark matter and dark energy implies that at the mo- 
ment we can not exclude the possibility of interactions 
in the dark sector. The confirmation of such a possibility 
would be of paramount importance because interactions 
between dark matter and dark energy could provide pos- 
sible solutions to the cosmological coincidence problem. 
In general, several papers have been published in this 
area (eg., Zimdahl, Pavon, Chimento 2001; Amendola et 
al. 2003; Cai & Wang 2005; Binder & Kremer 2006; Das, 
Corasaniti, & Khoury 2006; Olivares, Atrio-Barandela 
& Pavon 2008 and references therein) proposing that 
the dark energy and dark matter could be coupled. 

The aim of the present work is to investigate the 
observational consequences of the overall dynamics by 
using the A(t) cosmological model. Due to the absence 
of a physically well-motivated functional form for the 
A(t) parameter, we consider a power series form in H 
up to a second order. Doing so, we include the effects 
of the de-Sitter spacetime. The plan of the paper is as 
follows. The basic theoretical elements of the problem 
are presented in section 2 by solving analytically [for 
a spatially flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) geometry] the basic cosmological equations. In 
section 3 we place constraints on the main parameters of 
our model by performing a joint likelihood analysis uti- 
lizing the Union08 SNIa data (Kowalski et al. 2008), 
the shift parameter of the Cosmic Microwave Back- 
ground (Komatsu et al. 2009) and the observed Bary- 
onic Acoustic Oscillations (BAOs; Eisenstein et al. 2005; 
Padmanabhan, et al. 2007). Section 4 outlines the com- 
parison between the time varying vacuum model with 
the traditional A cosmology. Also, in section 4 we solve 
analytically the time evolution equation of the mass 
density contrast for the A(t) model while in section 5 
we present theoretical predictions regarding the forma- 
tion of the galaxy clusters. In section 6 we draw our 
conclusions. Finally, in the appendix we have treated 
analytically, the basic cosmological equations consider- 
ing that the time varying A(t) parameter can be ex- 
pressed with the aid of a power series expansion in H 
up to a third order. Note, that throughout the paper 
we use Ho = 70.5Km/sec/Mpc (Freedman et al. 2001; 
Komatsu et al. 2009). 



2 COSMOLOGY WITH A TIME 
DEPENDENT VACUUM 

In the framework, of a spatially flat Friedmann- 
Lemaitre-Robertson- Walker (FLRW) geometry the ba- 
sic equations which governs the global dynamics of the 
universe are 



Pm + PA 

and 

d(p m + pa) 
dt 



3H Z 



+ 3H(p m + P m + Pa + Pa) — 



(1) 



(2) 



where p m and pA are the matter density and vacuum 
density respectively, while P m = and Pa is the corre- 
sponding vacuum pressure. Note, that for simplicity we 
use geometrical units (8ttG = c = 1) in which pA = A. 



In order to study the above system of differential equa- 
tions we need to define explicitly the functional form 
of the A component. Within the framework of the A(t) 
model it is interesting to note that the equation of state 
takes the usual form of Pa = —pA(t) = — A(t) [Ozer & 
Taha 1987; Peebles & Ratra 1988]. 

On the other hand, introducing in the global dy- 
namics the above idea in a form of the time-dependent 
vacuum, it is possible to explain the physical properties 
of the dark energy. Considering now eq.(2), we have the 
following useful formula (see also Carneiro et al. 2008): 



p m + 3Hp„ 



-A 



and indeed, using eq.(l), we obtain: 
2H + 3H' 2 = A 



/ H dy p du t 

= J ~2 = 2 



(3) 



(4) 



(5) 



where the over-dot denotes derivatives with respect to 
time. Of course, the traditional A cosmology can be de- 
scribed by the above integration (eq.5) using a constant 
vacuum term A = const (for more details see section 
3.5). 

Now, from eq.(3), it becomes evident that in this 
cosmological scenario there is a coupling between the 
time-dependent vacuum and matter. Actually, the idea 
for possible interactions in the dark sector is not really 
new in this kind of studies. It has been shown that the 
coupling between dark matter and dark energy could 
provide possible solutions to the cosmological coinci- 
dence problem (eg., Zimdahl, Pavon, Chimento 2001; 
Amendola et al. 2003; Cai & Wang 2005; Binder & Kre- 
mer 2006; Das, Corasaniti, & Khoury 2006; Olivares, 
Atrio-Barandela & Pavon 2008 and references therein). 
In this context, one of the most important issues and 
unknowns is the precise functional form of the equation 
of state parameter w(a) where (a is the scale factor). 
The usual procedure is to derive a w(a) approximate 
functional form, by using a Taylor expansion around 
the present epoch (eg. Chevalier & Polarski, 2001; Lin- 
der 2003), which then provides approximate solutions of 
the global density evolution. However, the current ap- 
proach is somewhat different in the sense that we do 
not "design" the equation of state parameter such that 
to produce the desired (accelerated) cosmic evolution. 
Rather, we investigate whether a generalized vacuum 
component with w(t) = —1 and Pa = —ph{t) [Ozer & 
Taha 1987; Peebles & Ratra 1988] in the expanding Uni- 
verse allows for a late accelerated phase of the Universe 
and under which circumstances such a solution provides 
a viable alternative to the dark energy. 

Although, we do not have a fundamental theory 
to model the time-dependent A(t) function, we can 
parametrize the latter using a phenomenological ap- 
proach. Indeed, in a series of recent papers, authors (see 
for example Ray, Mukhopadhyay & Meng 2007; Sil & 
Som 2008 and references therein) have investigated the 
global dynamical properties of the universe considering 
that the vacuum energy density decreases linearly either 
with the energy density or the square Hubble paramo- 
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tcr. Also, Wang & Meng (2005), based on thermody- 
namics found that the vacuum energy density possibly 
decays as a power law. Alternatively, Carneiro et al. 
(2008) proposed a different pattern in which the vac- 
uum term is proportional with the Hubble parameter, 
A(a) oc H(a). However, this model fails to fit the cur- 
rent CMB data (see also section 3.4). In this context, at- 
tempts to provide a theoretical explanation for the A(t) 
have been presented in the literature (see Grande et al. 
2006 and references therein). These authors found that 
a time dependent vacuum could arise from the renor- 
malization group (RG) in quantum field theory. The 
corresponding solution for a running A(t) is found to 
be A(t) ~ ciH 2 (t) [where ci is a constant] and it can 
mimic the quintessence or phantom behavior and tran- 
sit smoothly between the two. It is worth noting that 
at late enough times the above solution asymptotically 
reaches the de-Sitter regime A ~ H 2 , as far as the global 
dynamics is concerned. 

In this paper, we parametrize the functional form 
of A(t) by applying a power series expansion in H up 
to the second order (see the appendix for a third order 
expansion which interestingly predict models with late 
acceleration but without initial singularities): 



A(t) =mH + n 2 H 2 



(6) 



Obviously, eq.(6) can be seen as a combination of the 
of the above ansatzs namely H(t) (Carneiro et al. 2008) 
and H 2 (t) [quantum field theory; Grande et al. 2006] re- 
spectively. It is now routine to integrate eq.(5) and ob- 
tain the Hubble function predicted by the current A(t) 
model: 

m e" lt/2 

= 7nF*7*rr ■ « 

where the range of /3(= 3 — 712) values for which the 
above integration is valid is (3 G (0, +00) (or n,2 < 3). 
Of course, if we consider different patterns for the vac- 
uum density then we can obtain different solutions for 
the Hubble parameter. Using now the definition of the 
Hubble parameter H = a /a, the scale factor of the uni- 
verse a(t), evolves with time as 



i(t) = a! (e nit/2 - l) 



2/0 



(8) 

where 01 is the constant of integration. As expected, at 
late enough times the above solution reduces to the de- 
Sitter universe. Note, that for /3 — > 3 and at early times 
the A(t) model tends to the Einstein de-Sitter case. Now 
from eqs.(7, 8) we can easily write the corresponding 
Hubble flow as a function of the scale factor 



H(a) 



ni 

13 



-13/2 



(9) 



Evaluating eq.(9) at the present time (a = 1) we obtain 

f3H 



Til 



1 + af 2 



(10) 



where Ho is the Hubble constant. From eqs.(9, 10), using 
the usual unit-less f2-parameterization, we have after 
some algebra that: 

H(a) 



E{a) 



{l-Q, m + Q, m a- p/2 ) 



(11) 
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Figure 1. Likelihood contours in the (fi m ,/3) plane. The 
contours are plotted where — 21n(£/£ max ) is equal to 2.32, 
6.16 and 11.83, corresponding to Icr, 2cr and 3cr confidence 
level. In Fig. la we present the likelihood contours that cor- 
respond to the SNIa, CMB and BAOs observational data. 
Finally, Figs, lb, lc and Id show the statistical results for 
different pairs. 



while the corresponding matter density parameter is: 
^m(cs) = f2 m a~ ,3 ' /2 /E(a). The normalized scale factor 
of the universe becomes 



a(t) 



t(o) = 



2/(3 



^ e f3(l-n m )H t/2 



1] 



2/(3 



/3(i - o, 

0/2 



-H^ln 



,/3/2 



E(a) 



(12) 



(13) 



where a\ l * = 57 m /(l — n m ). It is interesting to point 
here that the current age of the universe [a — 1, E(l) = 
1] is 



2 j \nQ„ 

t0 = p H ° 



i 



(14) 



We now investigate the circumstances under which an 
inflection point exists and therefore have an accelera- 
tion phase of the scale factor. This crucial period in the 
cosmic history corresponds to d(ti) = 0. Differentiating 
twice eq.(12), we then have: 

2/0 2 8 



a i 



W - 2)0, 



2(1 -On) 



/3(1 



which implies that the condition for which an inflection 
point is present in the evolution of the scale factor is 
P > 2. 



3 COSMOLOGICAL CONSTRAINTS 

3.1 The likelihood from the CMB shift 
parameter 

A very accurate and deep geometrical probe of dark en- 
ergy is the angular scale of the sound horizon at the 
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last scattering surface as encoded in the location of the 
first peak of the Cosmic Microwave Background (CMB) 
temperature perturbation spectrum. This probe is de- 
scribed by the so called CMB shift parameter (cf. Bond, 
Efstathiou & Tegmark 1997; Trotta 2004; Nesseris & 
Perivolaropoulos 2007) which is a normalized quantity 
and it is defined as 



R - 




dz 



(16) 



One of the merits of using the shift parameter in cos- 
mological studies is that its dependence on the Hub- 
ble constant is negligible (for details see Melchiorri et 
al. 2003; Nesseris & Perivolaropoulos 2007 and refer- 
ences therein). The shift parameter measured from the 
WMAP 5-years data (Komatsu et al. 2009) is R = 
1.71 ± 0.019 at zi s = 1090 [or a Ls = (1 + z^y 1 ~ 
9.17 x 10~ 4 ] and E(z) = H(z)/H is the normalized 
Hubble flow. Therefore, the corresponding Xcmb func- 
tion is simply written 



Xcmb(p) 



[R(p) - i-7i] s 

0.019 2 



(17) 



where p is a vector containing the cosmological param- 
eters that we want to fit. Note, that we sample the un- 
known parameters as follows: Q, m G [0.1, 1] and j3 £ [2, 5] 
in steps of 0.01. In Fig. la we present the la, 2a and 3a 
confidence levels in the (Q m , P) plane. It is evident that 
the f3 parameter is tightly constrained (/3 ~ 3.58) while 
the matter density parameter is not and all the values 
in the interval 0.1 < fi m < 1 are acceptable (see Ta- 
ble 1). However, following the WMAP 5-years results 
(Komatsu et al. 2009) of the full temperature perturba- 
tion spectrum AT/T, we can use an additional constrain 
which is n m h 2 = 0.1326 ± 0.0063. Thus, for h ~ 0.71 
(Freedman et al. 2001; Komatsu et al. 2009) we find 
0.24 < Q m < 0.29 (2a limits). 



constrain Q, m . In this case, the likelihood function can 
be written as: 

307 r 



XSNIafe) = ^ 



M th Kp)-M° bs (a0 



a, 



(18) 



where a; = (1 + z;) -1 is the observed scale factor of the 
universe, 2; is the observed redshift, n is the distance 
modulus fi — m — M — 51ogdL + 25 and dh(a, p) is the 
luminosity distance 



d L (a,p) = / A J 



(19) 



where c is the speed of light (= 1 here). Figure la also 
shows the la, 2a and 3a confidence levels in the (f2 m , /3) 
plane. Although, the (3 parameter is not constrained by 
this analysis the matter density parameter has an upper 
limit of f2 m < 0.29 within the la uncertainty (see Table 

!)■ 



3.3 The likelihood from BAOs 

In this section, we utilize the so called Baryonic Acoustic 
Oscillations. BAOs are produced by pressure (acoustic) 
waves in the photon-baryon plasma in the early uni- 
verse, generated by dark matter overdensities. Evidence 
of this excess has been found in the clustering proper- 
ties of the luminous SDSS red-galaxies (Eisenstein et 
al. 2005; Padmanabhan, et al. 2007) and it can provide 
a "standard ruler" with which we can put constraints 
on the cosmological models. In particular, we use the 
following estimator: 



A(p) = 



[zfB(a.)]V3 



/' 



da 



a 2 E(a) 



2/3 



(20) 



measured from the SDSS data to be A = 0.469 ± 0.017, 
where z a = 0.35 [or a s — (l + z s ) _1 ~ 0.75]. In this case, 
the Xbao function is given 

[A{p) - 0.469] 2 



Xbao(p) = 



0.017 2 



(21) 



3.2 The likelihood from the SNIa 

We now use the publicly available Union08 sample of 
307 supernovae of Kowalski et al. (2008) in order to 



* http://supernova.lbl.gov/Union. This catalog includes the 
following components in the error budget of the distance 
moduli: (a) <T to t, ( see Kowalski et al. 2008) obtained due 
to lensing, Milky way dust extinction and host galaxy pe- 
culiar velocities (b) the systematic error o- sys and (c) the 
uncertainty which is related with the light-curve fitting. For 
the latter uncertainty, Kowalski et al. take into account the 
stretch s and color c corrections via /is = mg— M+a(s— 1) — 
be for a specific cosmological model (Q m , w) = (0.29, —0.97). 
We would like to caution the reader that we do not min- 
imize x 2 over the parameters a and b, which implies that 
in the case of the A(i) model we may not be treating the 
third component of the error budget properly. However, wc 
would like to stress that according to Kowalski ct al. (2008) 
the corresponding constants a and b are rather insensitive to 
the assumed cosmological parameters (see their section 5.1). 
Thus, had we included the proper light-curve uncertainty in 
our fit wc would have obtained a larger solution space (see 
figure la). 
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Figure 2. The variance Ax 2 = X 2 — X^in aroun d the best 
fit Q m value when we marginalize over (i = 3 [A(t) oc H(t)]. 
The solid and the dashed line corresponds to the CMB and 
SNIa/BAO likelihoods respectively. 

It is evident (see figure la), that the matter density 
parameter is constrained (fl m ~ 0.28) by this analysis, 
while the j3 parameter is not (see also Table 1). 



to be more than 3a away, fl m — 0.80 (see solid line 
in figure 2), from the SNIa/BAO solution Q m ~ 0.32 
(see dashed line in figure 2). This implies that the func- 
tional form E(a) = 1 — Q m + f2 m a~ 3,/2 fails to fit the 
CMB data. We thus argue that the A(t) oc H{t) re- 
lation produces a discrepancy between the SNIa/BAO 
and CMB shift parameter which may lead to mislead- 
ing cosmological results. We further confirm the lat- 
ter result, by using a Bayesian statistics (see for ex- 
ample Davis et al. 2007), in which the corresponding 
estimator is defined as: BIC = \ 2 + klnN (where k 
is the number of parameters and N is the number of 
data points used in the fit). The next step is to es- 
timate the relative deviation between the two models 
ABIC = BIC niH - BIC niH+n2H \ In general a dif- 
ference in BIC of ABIC > 6, is considered strong evi- 
dence against that model which ocuurs the larger BIC. 
In our case, we find ABIC ~ 69 which implies a strong 
evidence against the A(t) oc H(t) model. 



3.4 The joint likelihoods 

We can combine the above probes by using a joint like- 
lihood analysis: 

Ctot{v) = £bao X £ cm b X £sNIa 



Xtot(p) 



2 2 2 

XBAO + Xcmb + ASNIa 



in order to put even further constraints on the param- 
eter space used. Note, that we define the likelihood 

estimator ^ as: Cj oc exp[— Xj/2]- The resulting best 
fit parameters, are presented in the second two rows 
of Table 1. The overall likelihood function peaks at 



0.29! 



P = 3.44 ± 0.02 and the correspond- 



ing Xtoti^m, P) is 310.2 (dof = 307). In this cosmolog- 
ical scenario the current age of the universe is found 
to be to ^ 14.1Gyr and the inflection point is located 
at {ai,ti) ~ (0.49, 0.44io )■ In Figs, lb, lc and Id we 
present, for various observational pairs, the correspond- 
ing likelihood contours (see also Table 1). 

Finally, it is worth noting that Carneiro et 
al. (2008) considered a different ansatz in order to 
parametrize the time dependence of the vacuum en- 
ergy. Their assumption is based on the fact that A(£) 
is proportional to the Hubble parameter (in our for- 
mulation ri2 =0). They found that this model fits the 
observational data (BAO+SNIa+CMB) at 2a level for 
Q, m ^ 0.43. In our case, if we marginalize over (3 = 3 
(or 712 =0), then the joint likelihood analysis provides 
a best fit value of f2 m ~ 0.35, but the fit is rather poor 
Xtoti^m) — 383 (dof = 308). We investigate a bit fur- 
ther this result and we reveal that the poor joint fit 
is due to the fact that the best fit value provided by 
the likelihood analysis of CMB shift parameter is found 

t Likelihoods are normalized to their maximum values. 
Note that the errors of the fitted parameters represent la 
uncertainties. 



3.5 The standard A-Cosmology 

In this section, we wish to remind the reader of some 
basic elements of the concordance A-cosmology in order 
to appreciate the differences with the A(t) cosmology. In 
the case of A = canst, it is straightforward to integrate 
eq.(5). Therefore, the Hubble function predicted by the 
A model is 



H(t) = J lea* § J±t 



(22) 



where A = 3H( 2 (1 — fi m ). Then the normalized Hubble 
function is written as 



E A (a) 



Hja) 
H 



[l-fi m + fi m a- 3 ] 1/2 



(23) 



while O m (a) = Q, m a 3 /E\(a). To this end, the scale 
factor of the universe is given by 



w = li-oJ smh 



OA 



tA[a) = ^7r=^^ lln 



VI - 0™ + E A (a) 



(24) 



(25) 



Comparing the A model with the observational 
data we find that the best fit value is Q. m = 0.28 ± 
0.02 with xL(^m) ^ 308.5 (dof = 308) in a good 
agreement with the 5 years WMAP data (Komatsu 
et al. 2009). Note, that Davis et al. (2007) using the 
Essence-SNIa+BAO+CMB and a Bayesian statistics 
found Q, m = 0.27±0.04, while Kowalski et al. (2008) uti- 
lizing the Union08-SNIa+BAO+CMB obtained Q. m = 

0.274±g:SEito.S" (for w 1). Obviously, our results 

coincide within la errors. 

The current age of the universe is given by 



3yi - a 



-.Ho 1 !*. 



VI - On + 1 



(26) 
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Table 1. Results of the likelihood function analysis. The I s * column indicates the data used (the last two rows corresponds to 
the inflection points). Errors of the fitted parameters represent la uncertainties. Finally, the current age of the universe to has 
units of Gyr (for H = 70.5Km/sec/Mpc). 
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SNIa 
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tr 



uncons. (Cl m = 
-,+0.12 
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3.44: 
3.44" 1 
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0.30 
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0.49 
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i /to 



0.29 
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0.33 
0.44 
0.44 
0.44 
0.44 



while the inflection point takes place at 



tlA 



3yi - a 



aiA = 



n/5 + 1 



1/3 



2(1 - a, 



(27) 



Therefore, we estimate toA — 13.9Gyr, tiA — 0.52toA 
and a/A ~ 0.58. Finally, using the previously de- 
scribed Bayesian statistics we find that ABIC = 
BIC n lH+n2 H 2 _ BI(J A ^ 5 Thig comparison i m pli es 

a preference for the usual A cosmology. 



4 COMPARISON BETWEEN DIFFERENT 
TYPES OF VACUUM 

In this section, we investigate in more detail the cor- 
respondence of the A(t) model with the traditional A- 
cosmology (see previous sections) in order to show the 
extent to which they compare. 



4.1 Compare the cosmic evolution 

Knowing now the parameter space (fi m ,/?) we present 
the evolution of the A(t) scale factor seen in the up- 
per panel of figure 3 as the solid line. It can been seen 
that it closely resembles the corresponding scale factor 
of the A cosmology (dashed line) . We have checked the 
cosmic phases of the A(t) scenario against the concor- 
dance cosmology by utilizing the deceleration parame- 
ter, q(a) — —(1 — adlnH/da). The evolution of the de- 
celeration parameter is presented in the bottom panel 
of figure 3, while in the insert figure we plot the rela- 
tive deviation of the deceleration parameter, A(q — oa), 
between the two vacuum models. We find the following 
phases: (a) at early enough times a < 0.21 the deceler- 
ation parameters are both positive with q > oa, which 
means that the cosmic expansion in the A(t) model is 
more rapid decelerated than in the A case, (b) between 
0.21 < a < 0.49 the deceleration parameters remain 
positive but q < oa, (c) then for 0.49 < a < 0.58 the 
traditional A model remains in the decelerated regime 
oa > but the A(t) is starting to accelerate q < 
and (d) for 0.58 < a < 0.80 the deceleration parame- 
ters are both negative and as long as q < oa the A(t) 
model predicts a much more acceleration than in the A 




Figure 3. Upper Panel: Comparison of the scale factor pro- 
vided by the A(t) model (solid line) with the traditional A 
cosmology (dashed line). Bottom Panel: The evolution of the 
deceleration parameter. In the insert panel we present the 
relative deviation A(q — (/a) of the deceleration parameters. 



model (the opposite situation seems to hold prior to the 
present epoch 0.80 < a < 1). In a special case where 
A(q — oa) — [q — qa, either at a ~ 0.21 or a ~ 0.80] 
the two vacuum models predict exactly the same ex- 
pansion of the universe. From figure 3 it becomes clear 
that the A(t) model reaches a maximum deviation from 
the A cosmology prior to a ~ 0.1 (z ~ 9) and a ~ 0.45 
(z ~ 1.2). Therefore, in order to investigate whether the 
expansion of the observed universe follows one of the 
above possibilities, we need a robust extragalactic dis- 
tance indicator at redshifts z > 1.2. Finally, the decel- 
eration parameters at the present time are go ~ —0.50 
and qoA — —0.57. 

4.2 Compare the linear growth factor 

In the framework of a time varying vacuum, the corre- 
sponding time evolution equation for the mass density 
contrast, in a pressureless fluid is given by (Arcuri & 
Waga 1994; see also Borges et al. 2008): 
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Table 2. Cosmological data of the growth rate of clustering 
(see Nesseris & Perivolaropoulos 2008). The correspondence 
of the columns is as follows: redshift, observed growth rate 
and references. 



.to 



bs 



Rcfs. 



0.15 0.51 ±0.11 

0.35 0.70 ±0.18 

0.55 0.75 ±0.18 

1.40 0.90 ±0.24 

3.00 1.46 ±0.29 

D + (2H + Q)D - 



Verde et al. 2002; Hawkins et al. 2003 
Tegmark et al. 2006 
Ross et al. 2006 
da Angela 2006 
McDonald 2005 



2HQ-Q 



D = 



(28) 



where p m = 3H 2 — A (see eq.l) and Q(t) = — A/p m . 
It becomes clear, that the interacting vacuum energy 
affects the growth factor via the function Q(t). Obvi- 
ously, in the case of a constant A [Q(t) = 0], the above 
equation reduces to the usual time evolution equation 
for the mass density contrast (see Peebles 1993). In this 
context, the growing solution as a function of redshift is 
given by: 




D A (z) 



E\{x) 



Ax 



(29) 



We now proceed in an attempt to analytically solve 
cq.(28). To do so, we change variables from t to a new 
one following the transformation 

y = exp(mt/2) < y < 1 . (30) 

Doing so, eq.(28) can be written: 

P 2 y(y-1) 2 D" +2f3(y-l)(5y-P)D' -2(6-/3)(/3-2y)£> = 0(31) 

where prime denotes derivatives with respect to y. We 
find that eq.(31) has a decaying solution for (3 < 8 of the 
form D\(y) = (y — l)^ 6 )/' 3 . The second independent 
solution (growing mode) of eq.(31) can be found easily 
from the following expression: 



Figure 4. Upper Panel: The evolution of the growth fac- 
tor for different vacuum models. The lines correspond to 
A(t) (solid) and A (dashed). Bottom Panel: Comparison of 
the observed and theoretical evolution of the growth rate 
of clustering f(z). Note, that data can be found in Nesseris 
& Perivolaropoulos (2008). The best fit for the A(t) model 
shows a ~ 10% difference from the traditional A model, es- 
pecially at large redshifts. 



D(y) = D 1 (y) 



(32) 



Inserting eq.(10) and eq.(13) into eq.(30), the y variable 
is related with the scale factor as: 



y = 



n m ) + o 

r, 



(33) 



In the redshift regime [a = (1 + z) - ] the combination of 
the above two equations lead to the following growing 
mode: 



D(z) = C(n m )(l + zf- 



/3)/2 



f 

■J z 



(x + 1) 



E 2 (x) 



where 



c(n m ) = 



- - n 



(2/3-4)//3 



■(34) 



(35) 



In the upper panel of figure 4 we present the 
growth factor evolution, derived by integrating eq.(29) 
and eq.(34), for the two vacuum models. Note that the 
growth factors are normalized to unity at the present 
time. Despite the fact that the global cosmological be- 
havior of the A(t) vacuum model is in a good agreement 



with the usual A cosmology (as it seen in figure 3), the 
two vacuum cosmological models trace differently the 
evolution of the matter fluctuation field. In particular, 
close to the present epoch (z < 0.3) the A(i) growth 
factor reaches a plateau, which means that the matter 
fluctuations are effectively frozen. It is obvious that the 
growth factor in the A(t) model is much greater than 
that of the concordance A cosmology. Indeed, assuming 
that clusters of galaxies have formed prior to the epoch 
of z ~ 1.4 (o ~ 0.42), in which the most distant clus- 
ter has been found (Mullis et al. 2005; Stanford et al. 
2006), the deviation (1 — D/Da)%, of the growth factor 
D(a) for the A(t) scenario with respect to the A solu- 
tion Da{o) is —51% while prior to the inflection point 
(ai ~ 0.5) we find —43%. We conclude that the behavior 
of the growth factor is sensitive to the different types of 
vacuum with D(z) > Da(z) and it is expected that this 
difference will affect also the predictions related with 
the formation of the cosmic structures (see section 5). 



4.3 Compare the growth rate of clustering 

We further compare the two vacuum cosmological sce- 
narios by utilizing the well known indicator of cluster- 
ing, namely the growth rate f(a) = dhiD/dlna (Pee- 
bles 1993). The corresponding parametrization of the 
growth rate of clustering can be achieved by introduc- 
ing a growth index 7 (see Wang & Steinhardt 1998) 
defined by 



/(o) = n^(o) 



(36) 
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Figure 5. Theoretical predictions of the fractional rate of 
cluster formation as a function of redshift. Note, that the 
solid and dashed line corresponds to A(i) and usual A cosmo- 
logical model respectively (erg = 0.80). Finally, for erg = 0.95 
we get: (i) the open triangles for A(t) and (ii) open circles 
for A. 



In order to quantify the growth index we perform a stan- 
dard x 2 minimization procedure (described previously) 
between the measured growth rate of the 2dF and SDSS 
catalogs (see Table 2; Nesseris & Perivolaropoulos 2008) 
with those expected in our spatially flat cosmological 
models 



X 2 (7) 



E 



fobs(zi) — /model («», 7) 



(37) 



where <Ji is the observed growth rate uncertainty. In 
the bottom panel of figure 4, we present the measured 
fobs(z) (filled symbols) with the estimated growth rate 
function f(z) = fl^ n (z) for the considered cosmological 
models. Notice, that for the A(t) model (solid line) we 
use (Q m ,P) = (0.29,3.44) and for the A case (dashed 
line) we impose Q m — 0.28. Also, in the insert panel of 
figure 4 we plot the variation of Ax 2 = X 2 (l) ~ Xmm(7) 
around the best 7 fit. We find that the growth index 
is 7 = 0.501^2 (x 2 / dof = 1-14) for the A(i), which is 
somewhat less (but still within la errors) than the A 
growth index, <y A = 0.62±° (x7 dof = 0.75). 



5 THE FORMATION OF GALAXY 
CLUSTERS 

In this section we attempt to briefly investigate the clus- 
ter formation processes by generalizing the basic equa- 
tions which govern the behavior of the matter pertur- 
bations within the framework of a A(t) flat cosmology. 
Also we compare our predictions with those found for 
the traditional A cosmology. This can help us to un- 
derstand better the theoretical expectations of the A(£) 
cosmological scenario as well as the variants from the A 
model. 

The concept of estimating the fractional rate of 
cluster formation has been brought up by different 
authors (cf. Peebles 1984; Weinberg 1987; Martel & 
Wasserman 1990 Richstone, Loeb & Turner 1992). The 
above authors introduced a methodology which com- 
putes the rate at which mass joins virialized structures, 
which grow from small initial perturbations in the uni- 



verse. In particular, the basic tool is the so called Press- 
Schechter formalism which considers the fraction of 
mass in the universe contained in gravitationally bound 
structures (such as galaxy clusters) with matter fluctu- 
ations greater than a critical value 8 C . Assuming that 
the density contrast is normally distributed with zero 
mean and variance a 2 (5) we have: 



dP(S) 



2lYCT 



exp 



2a 2 (M,z) 



dS 



(38) 



where 5 C is the linearly extrapolated density threshold 
above which structures collapse, ie, 8 C = 1.686. Note, 
that it has been shown that <5 C depends only weakly on 
flm (Eke, Cole & Frenk 1996). In this kind of studies it is 
common to parametrize the rms mass fluctuation ampli- 
tude at 8 /i _1 Mpc which can be expressed as a function 
of redshift as a(M,z) = crg(z) = D(z)ag. The current 
cosmological models are normalized by the analysis of 
the WMAP 5 years data cr 8 = 0.80 (Komatsu et al. 
2009). The integration of eq.(38) provides the fraction 
of the universe, on some specific mass scale, that has 
already collapsed producing cosmic structures (galaxy 
clusters) at redshift z and is given by (see also Rich- 
stone et al. 1992): 



erf 



V2a 8 (z) 



(39) 



Obviously the above generic of form eq.(39) depends on 
the choice of the background cosmology. The next step is 
to normalize the probability to give the number of clus- 
ters which have already collapsed by the epoch z (cu- 
mulative distribution), divided by the number of clus- 
ters which have collapsed at the present epoch (z = 0), 
F(z) — P(z)/P(Q). In figure 5 we present in a logarith- 
mic scale the behavior of normalized cluster formation 
rate as a function of redshift for the two cosmological 
models. In particular, for the traditional A cosmology 
we find the known behavior in which galaxy clusters 
appear to be formed at high redshifts z ~ 2 (see for ex- 
ample Basilakos 2003 and references therein), while the 
same general picture seems to hold for the A(t) model. 
However, in the latter case we find the following differ- 
ences: (i) clusters appear to form earlier (z ~ 3.5) with 
respect to the A model and (ii) prior to 2 ~ 0.4 the 
cluster formation has terminated due to the fact that 
the matter fluctuation field effectively freezes (see sec- 
tion 3.4). It is worth noting that the different formation 
rates between the two vacuum models, is due to the fact 
that the evolution of the corresponding growth factors 
are different (see the upper panel of figure 4). Finally, 
for a higher as value (as = 0.95) the corresponding clus- 
ter formation rate moves to higher redshifts [see figure 
5: A(£)-open triangles and A-open points]. The opposite 
situation is true for as < 0.80. 



6 CONCLUSIONS 

In this paper we study analytically and numerically 
the large and small scale dynamical properties of the 
FLRW flat cosmologies in which the "vacuum" energy 
is a function of the cosmic time A(t). Assuming that the 
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vacuum component can be expressed as a power series 
A = n\H + (3 — f3)H 2 , we find that the time evolu- 
tion of the basic cosmological functions are described 
in terms of exponential functions which can accommo- 
date a late time accelerated expansion, equivalent to the 
standard A model. Performing, a joint likelihood anal- 
ysis using the current observational data (SNIa, CMB 
shift parameter and B AOs) , we put tight constraints on 
the main cosmological parameters of the A(t) model. In 
particular, we find £l m ~ 0.29, (3 ~ 3.44 and the age of 
the universe is to ~ 14.1Gyr (for h = H /W0 ~ 0.705). 
Also, we compare the A(t) scenario with the traditional 
A cosmology. We find that the behavior of the global 
expansion in the A(i) model compares well with that of 
the usual A cosmology. However, there are differences 
especially when we consider the small scale dynamics. 
Indeed, we reveal that the A(i) cosmological model has 
two important differences over the considered A cosmol- 
ogy: 

• The amplitude and the shape of the linear growth 
of perturbations are different with respect to the A so- 
lution. As an example, prior to the inflection point the 
A(i) growth factor increases by a factor of ~ 43%. In 
this context, the growth index of clustering (7 ~ 0.50) 
is somewhat different with that of the A model (7A — 
0.62). 

• The large scale structures (such as galaxy clusters) 
form earlier (z ~ 3.5) with respect to those produced in 
the framework of the concordance A model (2 ~ 2). 



APPENDIX 

In this appendix we treat analytically, as much as pos- 
sible, the problem of the time varying A(i) parameter 
with the aid of a power series in H up to a third order: 
A = niH + n 2 H 2 + n 3 H 3 (n 3 / 0). The time evolution 
equation for the Hubble flow is obtained by eq.(5) as: 

H 

il n ? 

(40) 



dy 



y(n 3 y 2 - f5y + m) 



where /3 = 3 — 712 . In particular, the discriminant D — 
[3 2 — 4nin 3 characterizes the solutions of eq.(40) as: 

• Case 1: D > (f3 — An\ri2 > 0): The corresponding 
general solution of eq.(40) is written 



In 



Hem (H - Pi )pi(pi-p2) 



(H — pi) P2(P1-P2) 

pi, 2 
n2 = 3, pi 



nzt 
~2 



(41) 



where pi, 2 = - '2^ ^ 0- As an exam pl e f° r P — (or 
P2) we obtain 



Hit) 
and 



P2 



1 



a(t) = a\ 



1 + 



1 - \f\ - e" 3 "! 4 



(42) 



(43) 



where riz < and p 2 > 0. We would like to point out 
that as long as the cosmic time takes large values (t >• 



1), the A(t) model has the de-Sitter feature due to a(t) ~ 
e P2 '. On the other hand, it is very interesting the fact 
that this model has no initial singularity. Indeed, for 
t — > we get a(t) — > 01. 

Now, if we consider /3 7^ then the situation becomes 
complicated (see eq.41) but for the special case of pi = 
2p2 we can derive the following analytical solutions: 

P2 



H(t) = p 2 + 



and 



\J 1 - e" 3 "! 4 



a(t) 



oie' 



P2* 



1 + VI - e" 3P i* 



1 



p n 3 p 2 t 



(44) 



(45) 



where 713 < 0, p2 > and /3 < 0. Again, the A(t) model 
asymptotically reaches the de-Sitter regime a(t) ~ e 2p2t , 
while for t — > we again find no singularity a(t) — > 
ai. 

• Case 2: D — (f3 2 = 4nin 3 ): In this case the inte- 
gration of eq.(40) leads to the solution of: 



In 



H 



p 2 n 3 t 



H-p J H-p 



(46) 



where p = / 0. Now if f3 = (p = 0), which implies 
that rii = 0, then the solution of eq.(40) is given by 



— , n 3 < 

nzt 



(47) 



H{t)~- 
and 

a(t) = a-icv^ . 

• Case 3: D < {J3 2 - 4nm 2 < 0): In this case the 
integration of eq.(40) leads to the solution of: 



(48) 



In 



mH A 



+ 



2/3 



tan _1 G( J ff) 



n 3 H 2 -/3H + m J 
where G(H) = (2n 3 H - f3)/y/=D. 
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